
% find treated locations
treatedZ = GridsZ.tax > 1.15 ;

% employment DiD
sum( log( ( 1 - uZ ) .* popZ ./ ( ( 1 - uC ) .* popC ) ) .* treatedZ ) / sum( treatedZ )
    % 0.0430




